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Abstract 

For the two versions of the KdV equation on the positive half-line an initial-boundary 
value problem is well posed if one prescribes an initial condition plus either one bound- 
ary condition if q t and q xxx have the same sign (KdVI) or two boundary conditions if 
qt and q xxx have opposite sign (KdVII). Constructing the generalized Dirichlet to Neu- 
mann map for the above problems means characterizing the unknown boundary values in 
terms of the given initial and boundary conditions. For example, if {q(x, 0), q(0, t)} and 
{q(x, 0), q(Q, t), q x (0, t)} are given for the KdVI and KdVII equations, respectively, then 
one must construct the unknown boundary values {q x (0 7 t), q xx (0,t)} and {q xx (0,t)}, re- 
spectively. We show that this can be achieved without solving for q(x,t) by analysing a 
certain "global relation" which couples the given initial and boundary conditions with the 
unknown boundary values, as well as with the function fc), where satisifies the 

t-part of the associated Lax pair evaluated at x = 0. Indeed, by employing a Gelfand- 
Levitan-Marchenko triangular representation for the global relation can be solved 
explicitly for the unknown boundary values in terms of the given initial and boundary 
conditions and the function This yields the unknown boundary values in terms of a 
nonlinear Volterra integral equation. 

1 Introduction 

The Korteweg-de Vries (KdV) equation 

qt + cq x - Xq xxx + 6qq x = 0, c = 0, 1, A = ±1, (1.1) 
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appears in a wide range of physical applications. For example, it models one-dimensional 
small amplitude surface gravity waves propagating in a shallow channel of water. This 
equation has a celebrated history in the field of integrability; in particular it was for this 
equation that a method was first developed, later called the inverse scattering transform 
method, for solving the Cauchy (initial value) problem for integrable nonlinear evolution 
equations pQ, [2J- In 1968 Peter Lax realised that the key to the integrability of the KdV 
equation was the fact that it can be expressed as the compatibility condition of two linear 
eigenvalue equations now called a "Lax pair". Soon thereafter it was shown that a 
similar Lax pair formulation exists for the nonlinear Schrodinger equation 0|, as well as 
for the sine-Gordon and for the modified KdV equations. Actually, the existence of a 
Lax pair depending on a complex spectral parameter, denoted here by fc, is a defining 
property of integrability. 

The problem of solving the KdV equation on the half-line < x < oo, t > 0, was 
considered in using the general methodology introduced by one of the authors in and 
further developed in jjj. The starting point of this method is the simultaneous spectral 
analysis of the two eigenvalue equations that make up the associated Lax pair. This yields 
the solution of the KdV equation in terms of the solution of a matrix Riemann-Hilbert 
(RH) problem. This RH problem is formulated in the complex plane of the spectral 
parameter, the complex fc-plane, and its "jump" matrix is uniquely specified in terms of 
appropriate spectral functions, which are defined in terms of the initial condition q(x, 0) 
and the boundary values q(0, t), q x (0, t), q xx (0, t). Thus, the solution of a well posed initial- 
boundary value problem for the KdV equation can be expressed in terms of a matrix RH 
problem with known jumps, provided that one can construct the generalized Dirichlet to 
Neumann map, i.e. one can characterize the unknown boundary values in terms of the 
given initial and boundary conditions. It was shown in that this can be achieved, in 
principle, by analysing the so-called global relation, which is a simple algebraic relation 
coupling the spectral functions in the complex &;-plane. Since these functions are defined 
in terms of the initial condition and all the boundary values, the global relation provides 
an implicit characterization of the Dirichlet to Neumann map. A breakthrough in the 
analysis of the global relation was announced in [Hj , where it was shown that the global 
relation for the nonlinear Schrodinger equation can be solved explicitly in terms of the 
given initial and boundary conditions, as well as in terms of k), the solution of 

the i-part of the associated Lax pair evaluated at x = 0. Following this development, 
the generalized Dirichlet to Neumann map for several integrable nonlinear PDEs was 
analysed in [Hj, but the corresponding problem for equation (|1 . 1|> was not addressed. 
This omission was due to the fact that for the explicit solution of the global relation one 
uses an appropriate Gelfand-Levitan-Marchenko representation for k) and, in the 

case of equation jjl.ljl . such a representation was not known until now. In this paper, 
we present an appropriate Gelfand-Levitan-Marchenko representation for Hl.lfl and then 
use this representation to solve the global relation explicitly. 

The paper is organized as follows. In SectionElwe solve the global relation and present 
our main result. In Section we derive a Gelfand-Levitan-Marchenko representation for 
$W(t,Jfe). In Section H using the techniques developed in 0, we show that it is possible 
to solve the global relation for both versions of the KdV equation explicitly. Actually, by 
computing certain fc-integrals (see also |1 T| ^) we carry the technique of |S] one step further 
and show that the relevant expressions can be simplified substantially, see Section 
These results are discussed further in Sectional In Appendix A we review the general 
methodology of 0, namely we first analyse the "direct" spectral problem by defining 
appropriate spectral functions {a(k),b(k)} and {A(k), B(k)}, associated with the initial 
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condition and with the boundary values, respectively. We then formulate the "inverse" 
spectral problem through a matrix RH problem defined in terms of the above spectral 
functions. The solution of this problem yields the solution of the KdV equation with 
prescribed initial and boundary data, provided that the spectral functions satisfy the 
global relation. In Appendix B we discuss the generalized Dirichlet to Neumann map for 
the linear version of equation Hl.ljl . 

Remark 1.1. Equation is written in laboratory coordinates. We note that the 

choice A = 1 corresponds to the KdV equation with dominant surface tension. We also 
note that the usual caveat where one removes the q x term from equation 1)1.1(1 by chang- 
ing to travelling coordinates is not available in the quarter-plane without introducing 
a moving boundary problem (see |13|). If one drops the q x term arbitrarily (i.e. take 
c = in equation l|l.l[) ). the resulting initial-boundary value problem may be treated by 
a considerably simplified version of the analysis that is developed here. 

1.1 Notations and Assumptions 

We use the following notations and assumptions: 

• Subscripts with respect to x and t denote partial derivatives, for example 

dq d 3 q 

qt= di' qxxx = W 

• f(k) denotes the complex conjugate of the function f{k). 

• q{x, 0) = qo(x) denotes the initial condition. It is assumed that the initial condition 
belongs to the Schwarz space on the half-line, which will be denoted by 5(R + ). 

• We denote the boundary values at x — by 

q(0,t) = g Q (t), q x (0,t) = <&(*), q xx {0,t) = g 2 (t). (1.2) 

A subset of these functions is presribed as boundary conditions. We assume that 
these boundary conditions are sufficiently smooth and are compatible with qo(x) at 
x = t = 0. 

• Au,A 12 , A 2 i,A 22 denote the (11), (12), (21), (22) entries, respectively, of the 2 x 2 
matrix A. 

• D denotes the closure of the domain D in the complex fc-plane. 

• The dispersion relation Lu(k) is defined by 

oj{k) = 4Afc 3 + ck, A = ±l, c=0,l. (1.3) 

2 The Main Result 

Given the initial condition qo(x), define the spectral functions a(k) and b(k) by 

a{k) = V>2(0,A;), b(k) =<0i(O,fc), ImA: < 0, (2.1) 

where the vector (tpi(x, k), ^(x, k)) satisfies the following ODE, which is uniquely spec- 
ified in terms of q{x, 0): 

dipt , A , N/ , , x 
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where 

< x < oo, Imfc < 0, 

and 

lim ipi(x,k) = 0, lim ip 2 (x,k) = 1. (2.2b) 

3! — >OC X — >OC 

The global relation associated with the KdV equation is given by (see |S]) 

a(k)B(t,k) -b(k)A(t,k) = e 2luj{k)t c(t,k), Imfc < 0, < t < oo, (2.3) 

where w(fc) is defined by equation l|1.3|) . c(t, k) is an analytic function of k for Imfc < 
which is of 0(1 /k) as k — ► oo, and the functions {^4(i, fc), -B(t, fc)} are defined by 



i4(t,fe) = $ 2 (t,*), B(t,*) = -e 2M(t) '$i(t,fc), (2.4) 

where the vector ($i(i, fc), ^(i, ^)) T satisfies the following ODE, which is uniquely spec- 
ified in terms of {q(0,t),q x (0,t),q xx (Q,t)}: 

^ + 2tw(fc)*i = 2ikg (t)$ 2 + ffi(*)«2 + ^M*) - Ac ffo(^) ~ 2A .9o(i) 2 ) (*i + *a), 
^ = -2<fcso(t)$i +fl 1 (t)*i - ^ (#»(*) " AcgoW - 2Afl («) a )(*i + *2), 

(2.5a) 

where 

< t < oo, fc € C, (2.5b) 

and 

$i(0,fc) = 0, $ 2 (0,fc) = l. (2.5c) 

In the remainder of this section we give our main results and show that, for both 
versions of the KdV equation, the global relation H2.3JI can be solved explicitly for the 
unknown boundary values, ic. we give the generalized Dirichlct to Neumann map for equa- 
tion (| X . 1 f) . These results will be proved in Sections 03 and 0] by using a Gelfand-Levitan- 
Marchenko triangular representation for the spectral functions A(t, fc) and B(t, fc). 

2.1 The Gelfand-Levitan-Marchenko representation of $i and $ 2 

Proposition 2.1. The vector (t, fc), <E>2 (i, k)) T has the following Gelfand-Levitan- 
Marchenko representation 



VJ J- t \F2(t,a,k) 
where 

Fx(t,a,k) = ^-kB 1 + C 1 + ^so(t)Ai - i[£>! - ^9o(t)B 2 - g^Si(*Mi] 

F a (t, s, fc) = fc 2 A! + fcB 2 + C 2 + - [D t - ^rg {t)B 2 - _L ffl (t)Ai] , 

fc 4A oiA 
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and the six scalar functions {A\(t, s), B\(t, s), B 2 (t, s), C\(t, s), C 2 (t, s), D\(t, s)} satisfy 
the following well-posed characteristic-value problem: 

(d t + d s )Ax = 2ig () B 1 + 8iAL>i (2.7a) 
(fit - d s )B 1 = (G 2 - ^ 5o ) Ai - gi B 2 - 2ig Q C 2 (2.7b) 

(fit + d s )B 2 = (G 2 + ^5o)Ai - gxB x - 2ig C 1 (2.7c) 

(fit - d s )d = ~ (g + g ogi + gi) A x + (G 2 + ^9o) B i 

+ ^ 90 ~° 2 + 27a 5 ") B2 + 9l ° 2 (2Jd) 

(d t + 8 S )C 2 = -G 2 B X + (G 2 + 7^9 2 o)B 2 + gxC x + (2i + 2ig )D 1 (2.7e) 
(d t - d s )D x = -L( gi +gl- 2G G 2 + i s g)yl x + ^9od s B 2 

+ ^(<?o +gm)B 2 + (G a + ^9 2 o)Ci + G 2 C 2 , (2.7f) 
where go = 4fgo> <?i = 4f9i an d G 2 = %(g 2 — Xcgo — 2Ag§), with boundary conditions 
B x (t,t) = -ig (t), Cx(t,t) = i 5l (t), D x (t,t) = l -(g 2 (t) - Xcg Q (t) - 2\g (t) 2 ), (2.8) 



and 

Ax(t, -t) = B 2 (t, -t) = C 2 (t, -t) = 0. 
Proof. The proof of this Proposition will be given below in Section [3] 



(2.9) 



2.2 The generalized Dirichlet to Neumann map for KdVI 

For the KdV equation (| 1 . 1 1) with A = —1, an initial-boundary value problem is well 
posed if we specify the initial condition q(x,Q) = qo(x) together with a single boundary 
condition, say q(0,t) — go(t). In this case the analysis of the global relation yields two 
algebraic relations which can be solved explicitly to determine the unknown boundary 
values g\(t) and g 2 (t). We will show that 



OT 



1 



p+ ^±lA(t,p + )-p- 



n JdD P+ - P- 



(t, 2t - i)dr 

Hp-) 



' a(p+Y 



a(p-)' 



A(t,p-) dk, (2.10a) 



92(t) = -2g (t) 2 +tg (t)B 2 (t,t) - - ffl (t)Ai(t,t) 



3-2/3 



2 1 



+ - / J(k)e 
n JdD 



Ai'(C) iSidBx 
~2~dr~ 

Kp+) 

a(p+) 



o (*-r)V3 



2Bi 



(t,2T - t)dT 

, Hp-) 



a(p-)' 



A(t,p-) dk, (2.10b) 
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where A(t, k) is defined by the first equation in J23J, Ai{Q is a solution of Airy's equation 
given by 

Ai(Q = — [ e^ k+k3 /^dk, C = -3" 1/3 (i - r) 2 / 3 , (2.11) 
2tt Job 

the oriented contour dD is shown in Figure I2~T1 uj{k) is defined by equation (|1.3[) with 
A = — 1 and w'(k) denotes the derivative of w(fc), and p± are the two non-trivial roots of 
u>(p) — oj(k) 

P± {k)^- k -± l -{^-cY'\ (2.12) 



2.3 The generalized Dirichlet to Neumann map for KdVII 

For the KdV equation (|1.1J) with A = 1 , the initial-boundary value problem is well posed if 
we specify the initial condition q(x,0) = qo(x) together with two independent boundary 
conditions. To simplify the analysis, we consider the case where q(0,t) = go(t) and 
q x (0, t) = gi(t) are given and solve the global relation explicitly for the unknown boundary 
value (72 (i) ■ We will show that 

g 2 (t) = cg (t) + 2g {t) 2 - ig (t)B 2 (t,t) - ]- gi {t)A x {t,t) 



-1/3 



MO 



o (t-r)a/3 



o (*-r)V3 



G- 



Ai'(C) p3i dB x 
2~d7 

8iC\ -\ 



2B X 



(t,2r - t)dr 



dr 



3 

■9oW(o^T + 2 ^i) {t,2 T -t)dr 



' 2 di 



- { kio'{k)e- 2 ^ th -^-A(t,k)dk, (2.13) 
7T J dD a(k) 



where Ai(Q is given by (|2.11() but with C, — 3 _1 / 3 (i — t) 2 / 3 , uj(k) is defined by equation 
<|1.3|) with A = 1 and ui'(k) denotes the derivative of uj(k), and the contour dD, shown in 
Figure l2~2l is the oriented boundary of the domain 



L>i = {k G C, Im k < n Im w(k) > 0}. 
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Figure 2.2: The oriented boundary dD for the case A = 1. 



Remark 2.2. We note that the explicit solution of the global relation given in (|2.1()|1 
and (|2.13(l is expressed in terms of the given initial and boundary conditions and in terms 
of the kernel functions that appear in the Gelfand-Lcvitan Marchcnko representation of 
3?l(t, k) and <I>2 (£, k). We note that that it is also possible to express the Dirichlet to 
Neumann map in terms of $i(£, k) and & 2 (t, k) only (i.e. without reference to the kernel 
functions), although the resulting formulas are more complicated than those given here 
(see 0). 

Remark 2.3. If we take q(x,0) — as our initial condition then the spectral function 
b(k) defined in 12.1fl is zero and so the formulas given above in (|2.10() and (|2.13() simplify 
considerably (the fc-integrals in these formulas are now zero also). Furthermore, for the 
special case where we take c = in equation (| 1 - 1 1> - i.e. where the q x term is absent, the 
formulas H2.10|) and (|2.13(l simplify further in that, instead of the Airy function H2.11J) . 
they now involve the gamma function T(a) where a is independent of r. Thus, for 
instance, if q(x,0) = and we take c = and A = — 1 in equation then for KdVI 

the unknown boundary values gi(t) and g 2 (t) can be expressed in terms of q(0, t) — go(t) 
by the following formulas: 



9i(t) 
92(t) 



± go (t)A 1 (t,t) + -^r(l) 



1 



(i-r)V3 



dB l 



{t,2T-t)dr, (2.14a) 



-2.9oW 2 
1 



6l7T 



r(§) 



ig (t)B 2 (t,t) 
1 



(t - r) 2 /3 



-5i(t)Ai(t,t) 
3i dB l 



2 di 



2B, 



(t, 2r - t)dr, 



(2.14b) 



where {A 1 ,B ll B 2 , C 1 ,C 2 ,D 1 } satisfy the hyperbolic system lj!T7|l - (|2~g|) . 

Remark 2.4. It is shown in Appendix B that, in the linear limit of small qo(x), the 
expressions m l(2~Tn|l and f2~T3|l reduce to the same representations that are obtained by 
solving the global relation associated with the linearized KdV equation 



qt + cq x - Xq xxx = 0. 



(2.15) 
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3 The Gelfand— Levitan— Marchenko representation of 

In this section we verify Proposition 12 . II directly. Let the matrix fc) be defined in 

terms of ($i, $ 2 ) by 

*«(*,*)= (SEE *i(*.*)V (3.1) 

v ' ; ^i(i,fc) $ 2 (i,fc)y v ; 

Then equations (|2.5(l imply that this matrix satisfies the equation 



kQ (t) + Qi(i) + -Q 2 (t) 



$ (t) {t,k), (3.2a) 



where 

Qo(t) = G (t) (J , Qi(t) = gx(t) (I J) , Q 2 (t) = G a (i) , (3.2b) 

and 

G (t) = -2<a,(t), G 2 (i) = ~(g 2 {t) - Xcg (t) - 2\g (t) 2 ). (3.2c) 
Substituting the representation (cf. [H] and 0) 

$W(f,fc)=I + y (^k 2 A(t,s) + kB(t,s)+C(t,s) + ^D(t,s)^ e^W-^ds, (3.3) 

into equation Ij3.2a(l we obtain an integral equation involving the matrix-valued functions 
A, B, C , D. By multiplying by k 2 the following identity (obtained by integration by parts) 

k 3 F(s)e l ^ t - s ^ 3 ds = j (F(t)<r s - F{-t)a 3 e 2lu)(k)ta ^ 

-jj (kF{s) + iF s (s)a 3 )e iu ^ t - s ^ds, (3.4) 

we can eliminate the integral terms involving k . Then, using the identity obtained from 
multiplying l|3.4(l by k, as well as ()3.4|l itself, we can eliminate the integral terms involving 
fc 4 and fc 3 . Finally, by using the identity obtained from multiplying <|3.4[1 by 1/k we can 
eliminate the integral terms involving fc 2 . The remaining terms involve fc, 1, 1/k, and 
1/fc 2 . The latter term involves the product Q 2 D, which is equal to zero if the matrix D 
has the following form 

D(t,s) = D 1 (t,s)(~^ ~A. (3.5) 

Equating to zero the coefficients of the terms involving fc, 1, 1/fc we obtain the following 
equations: 

B t + a 3 B s a 3 = Q C + QiB + Q 2 A - -^Q A (3.6a) 

4A 

C t + a 3 C s a 3 - ^- (A t + a 3 A s a 3 ) = Q D + Q X C + Q 2 B 

- -L (Q B + Q X A) + -^QoA s a 3 (3.6b) 

b t + a 3 D s a 3 + — (A t ,<r 3 + <r 3 A ss ) = Q%D + Q 2 C + — (Q B S + Q x A s )a 3 , (3.6c) 
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Equating the coefficients of the boundary terms involving fc 2 , fc, 1, 1/fc we find the following 
boundary conditions: 

A(t, t) ~ a 3 A(t, t)a 3 = (3.7a) 

B(t,t) - a 3 B(t,t)a 3 = Q (3.7b) 

C(t, t) - a 3 C(t, t)a 3 =Q X - -^—Q A(t, t)a 3 (3.7c) 
D(t, t) - a 3 D(t, t)a 3 - jrriMt, *) CT 3 + ^A s (t, t)) 

= Q2- ^(Q B(t,t) + Q 1 i(t,t))a 3 ; (3.7d) 



and 



A(t, -t) + a 3 A(t, -t)a 3 = (3.8a) 

B(t, -t) + a 3 B(t, -t)a 3 = (3.8b) 

C(t, -t) + a 3 C(t, -t)cx 3 = ^Qoi(i, -t)a 3 (3.8c) 
-t) + a 3 D(t, -t)a 3 + ^(M*, -t)<r 8 + ^ s (i, -*)) 

= ^(QoB(«, -«) + Qii(i, (3.8d) 
We choose the matrices A,B,C in the following form, 

I(M) = ^(M)(J J), S(M)=(-^ 2 -^), C*(M)=(g g). (3.9) 

Substituting l|3.5|) and (|3.9|) into equations l|3.6f) we find the following differential equa- 
tions: 

(d t - d s )B 1 = (G 2 - ■^G )A 1 - 9l B 2 + G C 2 (3.10a) 

{d t + d s )B 2 = G 2 A 1 - 9l B x + GqCi (3.10b) 

(d t - d s )C x = ^-G d s A x - ■^- 9l A 1 + G 2 B 1 + (^-G - G 2 )B 2 
4i\ 4A 4A ' 

+ giC 2 -G D 1 (3.10c) 
(d t + d s )C 2 = ^(dt + d s )A 1 + (i^G - G 2 )B 1 + G 2 B 2 +g 1 C 1 - G Di (3.10d) 

(d t - d.)Di = ^idsAi - ~^G d s B 2 + G 2 {C 1 + C 2 ) - gi b x (3.10e) 

(d t + d s )Di = ^(dt + d s )d s A 1 + -^GodsB, + G^ + C 2 ) - 9l D x (3.10f) 

Furthermore, equations (|3.7|) and 13. 8|) yield the following boundary conditions: 

2B 1 (t,t) = G a (t) (3.11a) 

2C 1 (t,t) = 9l (t) - -LG (t)^i(i,i) (3.11b) 

2D 1 (t,t) = G 2 (t) - _L Sl (t)Ai(t,*) + ^-G Q (t)B 2 (t,t) (3.11c) 
4«A 4iA 

- -^(d t + d s )A 1 (t,t) = -G 2 + -±-G (i)£?i(M), (3.11d) 
4iA 4iA 
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and 

Ai(t,-t)=0 (3.12a) 

B 2 (t,-t)=0 (3.12b) 

C 2 (i,-t)=0 (3.12c) 

2D!(t, -t) - -^(d t + d s )A 1 (t, -t) = ^G Q B t {t, -t) (3.12d) 

= -^g 1 A 1 (t, -t) + ^G B 2 (t, -t). (3.12e) 

By subtracting (|3.10e() from (|3. lOff) . we find the equation 

{(9 t + d a )A x - giA x - G (B 1 + B 2 ) - 8iXD^ = 0. (3.10f) 

Similarly, by adding H3.11c(l and l|3.11d|l . we find the boundary condition 
(d t + d s )A 1 (t : t) - fli(t)Ai(t, t) + G (t)(fli(i, t) + B 2 (t, i)) - SiXD^t, t) = 0. (3. lid') 

In summary, we have derived the following result: Let Q^it, k) be expressed in terms 
of the matrices {A, B,C, D} by equation 1(3. 3|1 . Let these matrices have the form given 
by equations 1(3. 5 j) and Q3.9JI . Then k) satisfies equation H3.2J) provided that the 

scalar functions {A\, B\, B 2 , C\, C 2 , -Di} satisfy the differential equations l(3.10a|) -i (3.10e(l . 
<|3.10f I) together with the boundary conditions (|3.11a|) - (|3.11c|) . (|3. lld / l at s = t, and the 
boundary conditions (|3.12() at s = —t. 

We will now show that the above differential equations and boundary conditions are 
satisfied provided that the six functions {A±, B\, B 2 , C%, C 2 , Di} satisfy the differential 
equations l|3.1()a|l -l |3.1()e|l . the differential equation 

(d t + d.)Ai = - Gq{B\ + B 2 ) + 8iXD u (3.10F) 

as well as the boundary conditions ((3.11a() -l (3.11c() and (|3.12a|l - (|3.12c|) . Indeed, if equa- 
tion (|3.10f ')) is valid then equation ((3.1 Of \ is also valid. Furthermore, evaluating equation 
1 3.1 Of I) at s = t, it follows that equation 1(3.1 ld'f) is also valid. In addition, evaluating 
1 3.10f'() at s = —t and dividing the resulting equation by l/4iA we find 

2Cx(t, -t) -^{d t + d s )A 1 (t, -t) = ^-Go(B!(t, -t) + B 2 (t, -*)) - ^9iM(t, -*)> 

which, taking into consideration equations l(3.12a(l and l(3.12b(l . yields equation l(3.12d|l . 
Finally, equations (|3.12al) and (|3. 12b|) imply that equation (|3.12cll is identically satisfied. 

It is conve nient to replace (d t + d s )A 1 in the RHS of l(3.10d(l by the RHS of ((3.10f"(> : 
thus equation l(3.10d(l can be replaced by the equation 

(d t + d s )C 2 = ^ 9l A x - G 2 B X + (G 2 - ^-)B 2 + ^Ci + (2* - G )£>i- (3.10d') 
Equations l(3.11b(l and I3.11d'|l suggest the following change of variables 

C 1 (t,s) = C 1 (t,s)-^-G (t)A 1 (t,s) 
8iX 

Di(t,a) = Dx(t,s) - ^-gi^A^s) + ^-G (t)B 2 (t,s). 
8iX 8iX 

This yields the required result. 
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3.1 The spectral functions A(t, k), B(t, k) 

Equation (|3.1() . together with equation l|3.2(l and the integral representation for k) 
in equation l|3.3(l . imply the following expression for <f>i: 



$i(t,fc) = 



_ e -tw(k)t 



-fcS 1 (i,s)+c7 1 (i,s)-^ J D 1 (i,s) 



- 2fcBi(t, 2r - t) + 2Ci(t, 2r - t) + — So(*)Ai(t, 2r - t) 



Di(t, 2r - t) - -Lg (t)B 2 {t, 2r - t) - ^ ffl (t)A x (t, 2r - t) 



e 2ia,(*)(T-t) dTj 

e 2^(fc)(r-t) d 



where we have used the change of variables s i— » 2t — t. The definition of £?(i, fc) (equation 
lj2~H) ) then yields 

£?(t,fc) = 2fcB! - (2(?i + ^5o(*)ii) + ~(2£>! - j^9o(t)B 2 - -^ 9l (t)A x ), (3.13a) 



where 

N(t,k) = [ N(t,2r -t)e 2w{k)T dT. 
Jo 

An analogous expression can also be written for A(£, k). 



(3.13b) 



4 The generalized Dirichlet to Neumann map 

As was noted earlier (see also Appendix A), given qo(x) and a "proper" 1 subset of the 
boundary values {go (t) , g\ (t) , g 2 (t ) } we can characterize the unknown boundary values by 
the requirement that the spectral functions {A(t, k), B(t, k)} satisfy the global relation 
Ij2.3|) . For the Riemann-Hilbert problem formulated in Theorem A. 3 (see Appendix A 
below), the spectral functions {A(t, k), B(t, k)} are needed in the domain 

D 3 = {k S C,Imfc > Onlmw(fc) > 0}. 

However, since equation (I2.3|l is valid for Im k < we must transform the global relation 
from the lower half fc-plane into D 3 . We do this by using the invariance properties of the 
exponential e 2 M fe )*. 

For the case A = — 1, the two non-trivial roots of the dispersion relation uj{p) = ui(k) 
are such that k £ D3 => p + (k),p-(k) 6 Di where D\ is defined as 

Di = {k £ C,Imfc < Onlmcj(fc) > 0}. 

Evaluating the global relation (|2.3|l at p + and p^ we thus find two equations for the 
spectral functions A(t, k) and B(t, k), which are valid for k £ D3 

a{p+)B{t,p + ) - b(p+)A(t,p+) = e 2 ^ t c{t,p + ), (4.1a) 
a(p-)B(t,p-)-b(p-)A(t,p_) = e 2 ^ t c(t,p_), k £ D 3 . (4.1b) 



"Proper" here means those boundary conditions for which the initial-boundary value problem is at least 
linearly well posed. For example, if A = —1 the boundary value problem with q(0, t) given is well posed. 
Similarly, if A = 1 the boundary value problem with {q(0,t),q x (0,t)} given is well posed. 



IZ 
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These two equations give two independent relations coupling the known spectral func- 
tions {a(k),b(k)} with the unknown spectral functions {A(t,k),B(t,k)}. If we specify 
appropriate boundary conditions, for instance q(0,t) — go(t), then we can use equations 
(|4. 1|) to determine the unknown boundary values gi(t) and g-2 it) in terms of the given 
initial and boundary conditions. 

An analogous situation also exists for the case A = 1. In this case the two non-trivial 
roots of the dispersion relation uj{v) — uj(k) have the following properties: 



k E D 
ke D 



iM -^u+(k) 6flf, V-{k) G Di, 



where D ( 3 >, Dif> are shown in Figure 14.11 If we adopt the notation whereby v(k) — 
{v + (k), V-(k)} is defined such that k E D3 => v(k) g Di, then from (|2.3|l we have 

a(u)B(t, v) - b{v)A{t, v) = e 2iuj( - k)t c(t, u), k E D 3 . (4.2) 

This equation is a single relation coupling the known spectral functions {a(v),b(v)} with 
the unknown spectral functions {A(t, v), B(t, v)}. If we specify appropriate boundary 
conditions, for instance q(0,t) = go(t) and q x (0,t) = gi(t), then we can use equation 
1)4. 2|) to determine the unknown boundary value (72 (t) in terms of the given initial and 
boundary conditions. 

Before proceeding to solve equations (|4.1|l and l|4.2(l for the unknown boundary values, 
we first introduce certain integral identities that will be required for our analysis: 



OD 



<)D 



e 2Mk)(r-t') K ^ T ^ d7 



e 2iu ^ T -^K(t,T)d7 



dk = 



dk 



(3A)V3 



MO 

(f - T)V3 

Ai'(0 



2(3A) 2 / 3 J (f - r) 2 /3 



K(t,r)dT, (4.3a) 
K(t, r)dr, (4.3b) 



where Ai(Q is the solution of Airy's equation given in l)2.11[l . K(t, r) is a smooth function 
of the arguments indicated, 

< t' < t, 
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and the contour dD is the oriented contour shown in Figure l2~Tl In order to prove the first 
identity we first interchange the order of integration and then introduce new variables 

s = -2(3A) 1 / 3 (t' - r) 1 /^, C = (SX)- 1 / 3 ^' - r) 2 / 3 

to evaluate the k- integral. The result then follows from equation 1(2. 11(1 . Similar consid- 
erations apply to prove the second identity above. 

We will also make use of the following integral identities: 



i)D 



uo\k) 



w(fc) 



kui(k) 



no 



2iu{k){T-t') 



K(t,r)dT 



K(t,r)dT 



dk=TTK(t,t') 
i.ir 

dk - 



(4.4a) 



Ai(Q dK(t,r) 



dk 



2(3A)V3 J (f - t)!/3 Qt 



Ai'(C) dK{t : r) 



dr, 



(4.4b) 



4(3A) 2 /3 J (f - r) 2 /3 Qt 



dr. (4.4c) 



Identity (|4.4al) follows from the usual Fourier transform identity by mapping dD to the 
real axis. In order to derive expression 14.4b|) we first write the LHS as follows 



/ 


cu(k) f 


IdD 


Jo 



K(t,r)dT 



K(t,t') 



2/ 



dk. 



This ensures that the integral is well defined. We next perform integration-by-parts on 
the f-integral to obtain 



1 

dD 2i 



K(t,t') 



dk. 



The two terms involving K(t,t') cancel and, by interchanging the order of integration 
and using the identities in 1(4.3(1 . we can evaluate the fc-integral to obtain the expression 



l|4.4b(l . Similar considerations apply to prove the third identity 14.4c(l . 

4.1 The solution of the global relation for A = — 1 

We will now use the two global relations 1)4.1(1 to derive equations ((2.10(1 . i.e. an explicit 
expression for the unknown boundary values gi (t) and g 2 {t) . Substituting into equations 
(|4.1() the expression for B(t, k) given by l|3~T5|> and rearranging we find 



2C 1 --g Q A 1 =2( P+ +p-)B 1 - 



P+ -P- 



a(p^ 



a(p. 



1 



p+-p-" V a (p+) a (p~) 



and 



2D X + \g B 2 + = 2p+p^B 1 - , 

2 Ai p + -p- \a(p + ) 



b(P+) A(t, P+ )- b -^-A(t, P S 
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where {Ai, Z?2, Ci, D{\ are defined by equation l(3.13bjl . We first assume that a(k) ^ 
for k G D±. If we multiply these two equations by u)'(k)exp[— 2iu>(k)t'] and integrate 
along dD, the contour shown in Figure T2.ll then terms involving c(t,p±) will vanish 
since u>' (k)c(t, p±) / a(p±) is analytic and of O(l) for Imfc < 0, and the oscillatory term 
exp[2iuj(k)(t — t')] is bounded in D±. We use the fact that ku>'(k) = (3uj(k) — 2) together 
with the integral identities (|4.3(1 and 14.4(1 to obtain 



2nC 1 (t,2i/-t) = ng (t)Ai(t,2t'-t)-Tr3 



-1/3 



MO 



o (f-rjVs 



3^ ggi 
~2~dr 



-2Bi 



dD P+ - P- 



p + b ^A { t,p + )-pJ-HA {t ,p_) 



'a(p+Y 



a{p-)' 



(t,2T-t)dr 
dk, 



and 



2-nD x (t,2i! - t) 



~g {t)B 2 (t, 2t' -t) + ^gitfiAiit, 2t' - t) 



en 



B x (t,2t' - t) 



/' 



(t'-T)V3 



3dBi 
2~dV 



2iB 1 



(t, 2r - t)dr 



/{ky 



-2iuj(k)t' P+P- 
P+ -P- 



Kp+ \t, P+ )- b -^lA(t,p„) 



a(p+Y 



a(p-)' 



dk. 



Letting t' — ► t in these expressions and using the identities for Bi(t, t), Ci(i, t) and Di(t, t) 
in (ESI we find fZlfy . 

Finally, we note that if a(k) does have zeros for k € -Di, then the integral involving 
Lo'(k)c(t, k) exp[2iw(fc)(i — t')]/a(k) does not vanish, but rather it has a contribution due 
to the relevant residues. The simplest way to compute this contribution is to deform the 
contour dD to include the relevant zeros. 



4.2 The solution of the global relation for A = 1 

For the case of A = 1, we use the single global relation (|4.2j) to determine the unknown 
boundary value g2(t)- Substituting into equation l|4.2|l the expression for B{t, k) given by 
i|3.13|l and rearranging we find 

2A = \g,B 2 + - 2^ Bl + ^24 + -soil] + ^A(t, v) + e ^m c lfl, 

2 4i 1 2 J a(y) a{v) 

We first assume that a(v) 7^ for v € D\. If we multiply this equation by uj' (v)exp[— 2iu>(k)t'] 
and integrate along dDi, then the final term on the RHS vanishes because oj'(v)c(t, v)/a(v) 
is analytic and of 0(1) for Im k < 0, and the oscillatory term exp[2iui(k)(t — t')] is bounded 
in D\. Using the fact that ku)'(k) = (3cu(k) — 2) together with the integral identities (|4.3|) 
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and 14.4jl we find 



2TtDx{t, 2t' -t) = ^g (t)B 2 (t, 2t' -t) + ^g 1 {t)A 1 {t, 2t' - t) 



7T3- 1 / 3 



2 6 

Ai(0 



/' 



„ (f-r)V3 



(f - r)V3 



3i- 



dr 



3dB l 
2~dr 

4Ci+5d(t) 



(t, 2r - t)dr 



3i Mi 



At 



(t, 2t - t)dr 



OD 



a(k) 



Letting t' — > t in this equation and using the identity for Di(t, t) in (|2.8|1 we find l|2.13|l . 

As in the case for A = —1, we note that if a(k) does have zeros for k 6 D%, then the 
integral involving oj'(k)c(t, k) exp[2ia;(fc)(i — t')]/a(k) does not vanish. The simplest way 
to compute the contribution from this term is to deform the contour dD\ to include the 
relevant zeros. 



5 Conclusion 

The generalized Dirichlet to Neumann map for intcgrable nonlinear evolution PDEs can be 
characterized by the requirement that the spectral functions {a(k), b(k), A(t, k), B(t, k)} 
satisfy the so-called global relation. The Gelfand-Levitan-Marchenko representation of 
k) was used in |H] and [H] to solve the global relation explicitly for the nonlinear 
Schrodinger, the sine-Gordon, and the modified KdV equations (see also ^Hl)- We have 
shown in this paper that the global relation can also be solved explicitly for both cases 
A = ±1 of the KdV equation. 

The analysis of the KdV equation in the case A = — 1 is more complicated than 
the analysis of the case A = 1 since in the former case the single global relation Q2.3JI 
must be solved for two unknown boundary values. For economy of presentation we have 
concentrated on the case where <?o(^) is prescribed as a boundary condition and, by using 
certain invariant transformations of the dispersion relation uj{k) = — 4fc 3 + cfc, we have 
constructed explicit expressions for gi(t) and g2(t). The analysis is similar for the case 
that, say gi(t), is prescribed. 

The explicit solution of the global relation yields the unknown boundary values in 
terms of the given initial and boundary conditions and in terms of the kernel functions 
that appear in the Gelfand-Levitan-Marchenko representation of k). These kernel 

functions in turn satisfy a system of nonlinear ODEs which are equivalent to a system of 
nonlinear Voltcrra integral equations for which the nonlinearity now appears in a simple 
explicit form. The rigorous analysis of this system has not yet been carried out, however 
because the nonlinear integral equations are of Volterra type it should be straightforward 
to establish existence of the solution at least for small t or for "small" boundary conditions 

m 

We remark that the explicit characterization of the Dirichlet to Neumann map that we 
have presented here is a significant simplification on earlier attempts. Furthermore, this 
development appears to be particularly useful in the context of numerical simulations of 
soliton equations. Indeed, the result of [Hj has been used in ^2] for the implementation of 
the so-called method of nonreflecting boundary conditions for the numerical integration 
of the modified KdV equation. Actually, it was in J7| that the formulae of were 
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simplified by considering certain /c-integrals. In this paper we have extended this analysis 
and applied it to the KdV equation. 

Finally, we note that initial-boundary value problems for integrable nonlinear PDEs 
can also be analysed using PDE techniques |13 j -|15 | . In particular we note the remarkable 
result of [El where global well posedness is established. We consider our approach as 
complementary to the above approach: Having established a priori global existence for 
q(x,t) and therefore for {q(0,t),q x (0,t),q xx (0,t)}, our results can be extended globally. 
Therefore, our results imply an explicit characterization of the unknown boundary values 
which appears interesting both for analytical as well as computational considerations. 
Also, we note the Riemann-Hilbert formalism reviewed in Appendix A together with the 
powerful Deift-Zhou method ^Hj yield explicit results for the large i-asymptotics of the 
solution. 

Acknowledgement. This work was partially supported by the Department of Applied 
Mathematics and Theoretical Physics at the University of Cambridge, by Trinity College, 
Cambridge, and by the EPSRC. 

Appendix A Lax pair, spectral functions and the "in- 
verse" problem 

In this appendix we summarize the main results for the KdV equation presented in 0. 
Details and proofs can be found in [H] and 

Equation admits the Lax pair formulation 

fi x —ik[a s ,fi\ = Q(x, t, k)fi, (A. la) 

fit + iio{k)[a 3 , fi) = P{x, t, k)fi, k £ C, (A. lb) 

where uj(k) is defined by (|1.3|) , the cigenfunction fi(x, t, k) is a 2 x 2 matrix- valued function, 
and the functions Q{x, t, fc), P(x, t, k) are defined by 

Q(x,t,k) = ^(^2 -ia 3 )q(x,t), 

P(x,t, k) = -2ka 2 q{x,t) + <j x q x {x,t) + ^r(f2 - «cr 3 ) (q xx (x,t) - Xcq(x,t) - 2Xq 2 (x,t)) , 
with 




and [, ] denotes the usual matrix commutator 

[cr 3 ,A] = a 3 A - Acr 3 . 

A.l The "direct" spectral problem 

The spectral functions denoted by {a(k), b(k),A(t, k), B(t, k)}, are defined in terms of the 
initial condition and the boundary values by the direct map 

{q (x),g (t), 9l (t), g 2 (t)} ^ {a(k), b{k),A(t, k),B(t, k)}. 
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In order to set up the map qo(x) i— > {a(k),b(k)}, we first introduce the matrix-valued 
function ^^(x, k). This function is defined in terms of qo(x) as the following matrix 
solution of the x-part of the associated Lax pair (A.l) evaluated at t = 0: 

d x ¥ x \x,k) -ik[a 3 ,^ x \x,k)} = Q (x,k)^ x \x,k), (A.2a) 

where 

Qo(x,k) = —(a 2 - ia 3 )q (x), (A. 2b) 

lim y (x \x,k) = I, 0<x<oo, (A.2c) 

and for the first column of we take Im k > 0, whereas for the second column we take 
Imfc < 0. 

Definition A.l. The spectral functions a(k) and b(k) associated with the initial condition 
qo(x) are defined by 

a(k) = tp 2 {0,k), b(k) = V>i(0,fc), lmfc<0, (A.3) 

where ipi(x,k) and ip 2 (x,k) are the following components of the matrix-valued function 
~fy( x \x,k) defined by equation (A. 2), 



¥ x Xx,k)=(^ k ) (A4) 

ypi(x,k) ip 2 (x,k)J 

In order to set up the map {go(t), gi(t), g 2 {t)} >— > {A(t,k),B(t,k)}, we introduce the 
matrix-valued function <&W(i, k). This function is defined in terms of {go{t), g\{t), g 2 {t)} 
as the following solution of the i-part of the associated Lax pair (A.l) evaluated at x = 0: 

k) + kj(k)[* 3 , $ (t) (*, k)} = P {t, fc)* (t) (t, k), (A.5a) 

where 

P {t, k) = -2<j 2 g k + a igi + ^(cr 2 - ia 3 ) (g 2 (t) - Xcg (t) - 2Xg (t) 2 )^-, (A.5b) 

and 

$ ( * ) (0,fc)=7, (A.5c) 

and k G C. 2 

Definition A. 2. The spectral functions A(t, k) and B(t, k) associated with the boundary 
values {go(t), gi(t), g 2 (t)} are defined by 



' 2u f^B{t,k)\ _ /$i(t,jfc) 



A{t,k) J \<S> 2 (t,k)J> 



(A.6) 



where <&i(t, k) and & 2 (t, k) are the (12) and (22) components, respectively, of the matrix- 
valued function ^\t,k) defined by (A.3). 



2 If t — * oo, we take k such that Imu(fc) > for the first column of $ ( -*- ) , and k such that Imw(fe) < for 
the second column. 
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A. 2 The global relation 

Although the relationship between the initial condition and the boundary values of the 
solution of the KdV equation is complicated, this relationship takes a surprisingly simple 
form when expressed in terms of the corresponding spectral functions, see equation l|2.3[) . 
In the case < t < oo, equation (|2.3|1 becomes 

a(k)B(oo, k) - b(k)A(oo, k) = 0, Imfc < n Imw(fc) > 0. (A.7) 

The derivation of the global relation can be found in [5| . 



A. 3 The "inverse" spectral problem 

The inverse map 

{a(k), e- 2l ^ k)t b{k),A(T, k), e- 2l ^B{T, k)} ^ q(x, t), 
defines the function q(x, t) by 

q{x,t) = lim 4fc 2 [M(x,t,k)) 12 , < x < oo, < t < T, 

k — >oo 



(A.8) 



where M(x, t, k) satisfies the 2x2 matrix RH problem formulated in Theorem A. 3, which 
we state below without proof. For conciseness we give only the solitonless case. The 
inclusion of solitons is a simple procedure, the details of which can be found in jo] . An 
improved procedure for including solitons is presented in 

Theorem A. 3. j^J Let qo(x) 6 Define {a(k),b(k)} by equation (|A.3|I . Suppose 

that there exist smooth functions {go(t), 9i(t) 7 g2(t)} satisfying {gj(0) = d^qo^O)}^, such 
that the functions {A(t,k),B(t,k)}, defined by equation (|A.6|) . satisfy the global relation 



(|2.3|) . where c(t, k) is analytic in {k S C \ | Imfc < 0} and is of 0(l/k) as k — > oo. 3 
Let M(x, t, k) be a solution of the following 2x2 matrix RH problem: 

• M is sectionally meromorphic in k G C \ {k \ Imw(fc) = 0}. 

• Along the contours {k \ Imw(fc) = 0}, M satisfies the following jump conditions 



M_(.T,<,fc) = M + (x,t,k)J(x,t,k), lmuj(k) = 0, 



(A.9) 



where the matrices M±(x,t, k) are the limit values of M as k approaches {Imw(fc) = 
0}, and the jump matrix J(x, t, k) is defined in terms of the spectral functions a{k), 
b{k), A(T,k), and B(T,k) as follows: 



r(fc) 



2,1) 



J(x, t, k) = < 



-r(fc) e - 2i 
i 



i 



Vr(fc) 



-,2W 



keD 1 nD 2 
l-hKfc)| 2 

- 7 (fc)e 2ie 

kelhnD; 



7 (fc)e~ 2if 
1 



(A.10) 



-T{k)e- 2W 
1 



k e D 3 n D 4 , 



3 For the case < t < oo, equation <12.'M is replaced by ljA.7fl and we further assume that go(t), gi{t), g2(t) 
belong to 5(K+). 
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where 

0(x, t, k) = -kx + Lo(k)t; (A.ll) 
the functions j(k), T(k) are given by 

7 (A) = 2H=, fce^n^, (A.12a) 

a(k) 

r(fc) = -^4^=, keTh^Th, Thnth, (A.i2b) 

o(fc) d(fc) 



d(k) = a(k)A(T,k)-b(k)B(T,k), k e D 2 ; (A.12c) 
one? i/ie domains D\-D± are defined as: 



D 1 = {k G C, Im jfe < H Im > 0}, 
D 2 = {fe G C, Im jfe < n Im u(k) < 0}, 
£» 3 = {k G C, Im A; > n Im to(k) > 0}, 
Di = {k G C, Im A; > n Im < 0}. 



(A.13) 



, k 

• M(x, t, k) has a pole at k = satisfying 



M(x,t,k) = I + [- I , fc^oo. (A.14) 



i/(a:,t) /0 1 



Then we have the following: 

(1) M(x, t, k) exists and is unique. 

(2) q(x,t) defined in terms of M(x,t,k) by equation (|A.8J) satisfies the KdV equation 

CD- 

(3) q{x,i) satisfies the initial condition q(x,0) = qo(x), and furthermore, it has the 
following boundary values: 

q(fl,t) =9o(t), q x (0,t) =gi(t), q xx (0, t) = g 2 {t). 
Proof. Details of the proof of this theorem can be found in 0. 

Appendix B The linear limit 

In this appendix we give the generalized Dirichlct to Neumann map for the linear version 
of equation Hl.ljl . We also prove that, in the linear limit, our construction of the Dirichlct 
to Neumann map for the KdV equation in Section [3 reduces to the corresponding map 
in the linear case. 

In the approximation of small qo, go, gi, gi (or small go and small t) then, from 
Proposition ^. II we find Bi(t,s) = Bi(t + s) and thus Bi(t,t) = Bi(2t) = —igo(t). Hence 



B 1 {t,2T-t) = B 1 {2r) = -tflb(r). 
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Similar considerations apply for the other kernel functions and so we have 
B 1 (t,2T-t)^-ig {T), 

Cl(t,2T-t)w^l(T), 

Di(t,2r-t)« J[fla(r)-A<go(T)], 



Ax{t,2r 


-t) p 


3 


B 2 (t,2r 


-*)<= 


a 


C 2 (t,2T 


-i)p 


a 



and, from (|2.2p. 



o(fe) - 1, b(k) 



B.l The linear limit for KdVI (A = -1) 

Taking the linear limit of 1(2.10(1 and substituting the above expressions, we find 



3i(f)»- / iku'(k) \ / e 2i "^-Vg Q (T)dT 



dk 



1 



to' (k)e- 2w{k)t 

dD P+ - V- 



p+qo(p+) -P-Qoip-) 



dk, 



and 



g 2 (t) » ~ / fcV(fc) / e^W^^dr 

2i 



dk 



+ 



-2iu(k)t P+P- 



u'(k)e 
od P+ - P- 



q (p+) - qo(P-) 



dk, 



where 



Qo(k) 



e- 2lkx q (x)dx. 



By using the identities (|4. H|> and (|4.4(l we see that these are the same formulae we get by 
solving the global relation associated with equation 1(2.15(1 with A = — 1, see \12\ . 

B.2 The linear limit for KdVII (A = 1) 

Taking the linear limit of 1(2.13(1 and substituting the above expressions, we find 



92(t) « cg (t) + - / k 2 Lo\k)\ e 2 ^ k ^g (T)dr 



71" Jd 



dk 

kcu'(k) I e^WCr-t] < gi ( T )dT 

Ai 



dk 



- I u'ifye-^Wqoifydk, 

K J 3D 



where 



<?o(fc) = / e- 2ikx q (x)dx. 



By using the identities 14.3(1 and ((4.4(1 we see that these are the same formulae we get by 
solving the global relation associated with equation 1(2.15(1 with A = 1, see |12( . 
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